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ABSTRACT 


A fully  developed,  two  dimensional  turbulent  flow  is  established 
in  a water  channel  with  a free  surface.  The  unstratified  case  is  first 
examined  in  regards  to  requirements  for  two  dimensionality,  steady  state, 
and  properties  near  the  free  surface  with  hot  film  anemometry.  The  flow 
is  then  stratified  by  injecting  saturated  salt  water  at  mid  depth.  Down- 
stream from  the  salt  injector  in  the  test  section  the  properties  of  the 
flow  are  then  examined.  Cross  sections  in  the  center  of  the  channel  one 
foot  apart  are  taken  of  the  mean  and  streamwise  fluctuating  velocities 
as  well  as  the  mean  density  profiles.  This  information  along  with  the 
boundary  conditions  of  the  stress  at  the  wall  and  the  pressure  gradient 
at  the  wall  are  then  fed  into  the  numerical  program  of  the  first  integral 
of  the  equations  of  motion.  This  yields  the  momentum  and  salinity  fluxes 
as  well  as  the  inertial  terms  and  pressure  gradient  as  a function  of  depth 
and  the  local  Richardson  Number.  The  salinity  flux  and  Austausch  coeffi- 
cient behave  with  the  expected  reciprocal  relationship  to  the  Richardson 
Number.  The  momentum  results  do  not.  The  flow  exhibits  internal  regimes 
in  which  the  dynamics  are  different.  In  the  lower  region  of  the  flow 
the  turbulent  momentum  flux  is  negative,  u’  and  w'  are  either  both  pos- 
itive or  negative,  and  the  mean  gradient  steepens  in  the  streamwise 
direction.  Also  the  turbulent  momentum  flux  is  balanced  by  the  mean 
advective  fluxes.  The  possibility  of  the  critical  layer  mechanism  trans- 
ferring momentum  from  internal  waves  to  the  mean  flow  is  considered  in 
the  last  section.  In  the  upper  region  of  the  flow  the  turbulent  mo- 
mentum flux  is  positive  and  is  balanced  by  the  vertical  integral  of  the 
streamwise  density  gradient.  The  momentum  Austausch  coefficient  is 
large  as  is  the  local  Richardson  Number.  This  is  probably  due  to  the 
fact  that,  the  local  Richardson  Number  as  defined,  does  not  properly 
assess  the  turbulent  kinetic  energy  available  for  mixing.  Energy  is 
advected  into  the  region  as  the  flow  is  developing. 


1.  Introduction 


One  area  of  oceanography  that  has  attracted  considerable  attention 
is  the  existence  of  seasonal  thermoclines,  main  or  permanent  thermoclines, 
and  small  scale  microstructure  or  layering.  This  paper  is  concerned  with 
one  mechanism  that  might  be  responsible  for  these  phenomenon;  that  is, 
the  effect  upon  shear-induced  turbulence  of  a stable  density  gradient. 

A fully  developed,  two  dimensional  classical  turbulent  flow  is 
established  similar  to  that  of  Laufer  (1951)  and  is  then  stratified  by 
injecting  saturated  salt  water  into  it.  This  approach  was  suggested  by 
the  experiments  of  Ellison  and  Turner  (1960),  although  it  has  been  sig- 
nificantly extended  in  parameter  space.  Because  the  state  of  technology 
prevented  them  from  including  advection,  which  is  central  to  the  dynamics, 
their  small  bulk  Richardson  Numbers,  and  their  small  aspect  ratio  (their 
flow  is  not  two  dimensional),  their  results  are  preliminary,  and  the 
interesting  features  of  this  experiment  are  missed. 

In  the  test  section,  which  is  sufficiently  downstream  of  the  in- 
jecting device  to  avoid  Its  effect,  profiles  of  ^(z),  0(z)  and/l?(z) 
are  taken  at  cross  sections  one  foot  apart.  This  information  along  with 
the  boundary  conditions  of  wall  stress  and  the  pressure  gradient  at  the 
wall  are  fed  into  a numerical  computer  program  of  the  first  integral 
of  the  equations  of  motion.  From  this  the  momentum  and  mass  fluxes  as 
well  as  the  advective  and  pressure  gradient  terms  and  local  Ri.  No.  as 
a function  of  the  normalized  depth  are  obtained  for  various  overall 
Ri.  Nos.  The  conservation  of  mass  across  the  two  sections  is  also 
determined  as  a check  on  the  results. 

Unfortunately  in  attempting  to  model  oceanographic  mechanisms  in 
the  laboratory,  compromises  were  necessary.  The  ocean  currents  exhibit 
quite  large  Re.  Nos.  and  Ri.  Nos.  In  going  from  the  ocean  to  the  labor- 
atory, the  physical  scale  changes  by  several  orders  of  magnitude.  For 

a constant  Re.  No.,  the  Ri.  No.  goes  as  the  cube  of  this  physical  scale. 

For  this  reason  the  Re.  No.  was  limited  to^'lO.OOO  and  the  range  of  the 

bulk  Ri.  No.  was  up  to » 0.40. 


2.  Experimental  Apparatus  and  Procedures 
a)  The  Water  Channel 

A twenty-six  foot  plastic  flume,  eight  Inches  wide,  was  constructed, 
supported  by  wooden  frames.  The  object  was  to  create  In  the  neutral  case  a 
fully  developed,  two  dimensional,  steady  state  turbulent  flow,  the  analog  of 
Laufer's  wind  tunnel  experiment,  (1951),  which  he  thoroughly  Investigated.  In 
this  turbulent  system  the  turbulent  energy  goes  through  all  of  Its  stages  of 
transformation  across  the  channel  section  - turbulent  energy  production  from 
the  mean  flow,  energy  diffusion,  and  turbulent  and  laminar  dissipation.  From 
the  test  section  the  flow  proceeds  over  a dam  at  the  channel's  end.  This  con- 
trols the  depth  of  the  flow  for  a fixed  flow  rate  and  flow  rate  Re. No.  The 
fluid  then  recirculates  through  a 2000  gallon  sump  which  acts  as  a thermal  and 
salinity  sink.  A bronze  pump  with  a Reeves  variable  speed  pulley  then  Issues 
the  fluid  through  an  Induction  type  flow  meter.  The  packing  gland  on  the  pump 
Is  adjusted  so  that  no  air  may  enter.  From  the  flow  meter  the  fluid  passes 
through  a full  flow,  bronze,  five  micron  filter  and  thence  to  the  settling  tank 
at  the  channel's  mouth.  In  the  settling  tank  the  fluid  passes  through  a large 
mesh,  a large  honeycomb,  a contraction  section,  and  finally  through  a small 
honeycomb  to  the  channel's  beginning.  The  details  of  the  apparatus  are  shown 
in  figure  1.  The  flow  was  stratified  by  placing  five  salt  Injectors  fourteen 
feet  downstream  from  the  channel  entrance.  At  this  point,  the  ratio  of  distance 
from  the  entrance  to  depth  was  100.  The  flow  was  fully  developed  turbulence  -- 
the  finite  amplitude  question.  The  test  section  was  five  feet  downstream  from 
the  injectors  where  their  wake  had  dissipated  and  the  flow  had  returned  to  steady 
state  conditions  for  the  neutral  case.  The  test  section  was  six  feet  upstream 
from  the  dam  to  avoid  any  end  effects,  such  as  Internal  hydraulic  jumps. 
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The  saturated  salt  water  injection  rate  was  regulated  by  the  con- 
trol tank.  In  the  control  tank  a float  attached  to  a lever  operated  a mercury 
level  switch  which  opened  and  closed  a valve  from  a storage  tank  of  445  gal- 
lons capacity.  The  depth  of  water  in  the  control  tank  was  regulated  to 
+ 

0.25  inches.  The  constancy  of  the  injection  flow  rate  is  shown  in  figure 
3.  The  injectors  were  placed  at  a normalized  distance  from  the  bottom  of  the 
channel  of  xhe  injection  flow  rate  was  6cm. /sec.  and  the  jets  lost 

their  definition  after  8cms.  from  the  injectors. 

b)  Density  Measuring  Techniques 

The  density  profiles  were  measured  by  taking  samples  of  the  fluid 
at  various  depths  and  feeding  th«n  into  a precision  hydrometer.  The  density 
probe  had  a horizontal  dimension  of  0.2  inches,  and  a vertical  dimension  of 
0.035  inches,  and  traversed  the  flow  on  a vernier  carriage.  The  hydrometer 
was  a Brooks  Instrument  continuously  reading  type  where  the  liquid  flowed  in 
through  a needle  valve  at  the  bottom  and  out  through  two  discharge  tubes  in- 
side the  glass  cylinder.  The  hydrometer  was  calibrated  by  measuring  the  den- 
sity of  fresh  water  at  several  different  temperatures  and  comparing  the  re- 
sults with  the  data  from  the  Handbook  of  Chemistry  and  Physics.  The  resolu- 
tion was  - .0003  gm./cc.  This  instrument  is  shown  in  figures  4 and  a discus- 
sion of  its  use  is  presented  in  secllon  four. 

c)  Velocity  Measuring  Techniques 

The  measurement  of  the  velocity,  a fluctuating  vector,  in  a vary- 
ing density  field  is  a considerably  more  difficult  proposition,  and  a great 
deal  of  time  went  into  developing  this  instrumentation.  A pitot  tube  was 
used  in  the  conventional  manner  for  calibration  and  for  mean  velocity  profiles 
in  the  neutral  case.  It  proved  useless  in  the  stratified  case  as  figure  6 
indicates.  The  theory  of  a pitot  tube  requires  that  the  fluid  decelerate 
isentropically  to  the  stagnation  point.  Considering  a turbulent  flow,  we 
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integrate  the  momentum  equation  along  the  center  stream  line  from  the  point 
where  the  velocity  is  that  of  the  flow  to  the  stagnation  point  where  the 
velocity  is  zero  (about  four  diameters  of  the  tube,  Goldstein,  1965). 

1-)  ^ ~~  along  the  streamline 

Now  let: 


where 

and 


if 


where  U»  is  the  velocity  of  the  flow  at  that  point. 

Then,  equation  1.)  becomes: 

2.)  — I if  - ^ ^ u/W*  ^ ^ 

^ %u  ^ fp  j ^ 


Rearranging: 


3.)  4. 
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over  most  of  the  flow. 

then: 

4.) 
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5.) 
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The  pitot  tube  can  be  used  in  a neutral  turbulent  flow  for  mean  velocity 
profiles  and  calibration. 

In  as  much  as  we  want  to  measure  turbulent  velocities  in  a strati- 
fied field,  we  need  an  instrument  that  has  a high  resolution  in  space  and  time 
and  is  not  effected  by  a varying  density  field.  This  is  the  hot  film 
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anemometer.  Hot  wires  have  been  used  successfully  for  years  in  air  born 
turbulence.  In  going  from  air  to  water  born  turbulence  five  major  problems 
are  encountered  which  have  only  recently  been  solved.  Even  then  hot  films 
are  tempermental  in  a liquid,  particularly  a salty  one.  Dirt  in  the  water 
building  up  on  the  hot  film,  ground  loops  creating  leakage  from  the  hot  film 
through  the  conducting  fluid  to  the  channel  supports,  air  bubbles  coming  out 
of  solution  and  attaching  to  the  film,  and  deterioration  of  the  insulating 
quartz  film  and  hydrolizing  of  the  platinum  heating  element  are  initially 
encountered.  Also,  the  low  over  heat  ratio  necessary  to  minimize  the  air 
bubble  problem  causes  the  film  to  be  a good  thermometer.  Small  temperature 
changes  of  the  fluid  (^1°C)  cause  a large  change  in  the  calibration  curve. 

To  cope  with  these  problems  a full  flow,  five  micron  filter  was  installed  in 
the  recirculating  line  and  the  channel  was  covered  with  a polyethylene  sheet. 
This  completely  eliminated  the  dirt  problem.  The  channel  was  made  of  1/2  inch 
plastic  and  the  supports  were  wooden,  thus  eliminating  the  ground  loop  problem. 
The  2000  gallon  underground  sump  allowed  the  water  to  come  to  equilibrium  in 
regards  to  dissolved  air  and  minimized  the  air  bubble  problem.  With  the  above 
conditions  the  life  time  of  the  films  was  usually  about  seven  hours  in  salt 
water,  sufficient  time  to  complete  a run.  In  fresh  water  the  life  time  was 
much  longer.  The  films  exhibited  a shelf  life  after  having  been  used  in  salt 
water.  The  2000  gallon  sump  was  a thermal  sink  and  kept  the  temperature  of 
the  fluid  constant,  usually  to  better  than  0.1“C  during  the  running  time  of 
three  hours  for  the  velocity  measurements.  This  caused  no  drift  in  calibra- 
tion. 

The  hot  film  is  shown  in  figure  5 and  the  electronic* in  figure  2. 
The  hot  film  traversed  the  center  cross-section  on  a vernier  carriage.  In  the 
region  near  the  wall  where  there  are  strong  gradients  an  ocular  micrometer  was 
used  to  determine  the  distance  from  the  wall.  The  instrument  was  aimed  at  a 
shallow  angle  to  the  channel  bottom  and  the  distance  between  the  hot  film  and 
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its  reflection  in  the  bottom  was  measured.  The  micrometer  read  to  .0001 
inch  and  the  resolution  of  the  distance  was  .0005  inch.  The  hot  film 
and  its  electronics  are  Thermo-Systems  equipment  of  the  constant  tempera- 
ture variety,  reducing  the  non-linear  problem  in  the  region  of  large  fluc- 
tuations. For  a fluctuating  to  mean  velocity  ratio  of  .20  the  fluctuating 
velocity  error  is  only  a few  per  cent  {Hinze,  1959).  The  hot  film  and 
electronics  response  is  flat  out  to  5 KHz,  well  beyond  the  spectral  range  of 
water  born  turbulence.  The  hot  film  is  made  by  drawing  a quartz  rod  to  one, 
two  or  six  mils.  The  one  and  two  mil  rods  were  used  in  the  neutral  case; 
whereas,  the  six  mil  was  used  in  the  stratified  case  beacuse  of  its  better 
durability.  An  80  angstrom  platinum  covering  is  sputtered  on  the  center 
region  of  the  quartz  rod  and  gold  is  sputtered  on  the  ends  to  be  welded  to 
the  supporting  needles.  The  corner  weld  sections  and  the  gold  are  covered 
with  epoxy  resin  while  the  platinum  heating  surface  is  insulated  by  a few 
angstroms  of  sputtered  quartz.  Because  of  the  platinum  thickness  to  length 
ratio,  there  was  no  problem  of  end  losses.  The  film  lengths  varied  from  0.02 
inch  to  0.08  inch. 

d)  Calibration  of  the  Hot  Film 

The  heat  transfer  from  a hot  film  is  more  complicated  than  from 
a hot  wire.  It  is  a function  of  many  variables  some  of  them  not  easily  define- 
able.  The  several  heat  transfer  theories  from  heated,  infinite  cylinders  were 
of  little  help  as  the  results  did  not  collapse  to  a single  curve.  Because  of 
this,  empirical  calibration  was  necessary. 

Two  methods  were  employed.  First,  two  concentric  plastic 
cylinders,  electrically  insulated  from  ground,  were  filled  with  the  desired 
fluid  and  brought  up  to  solid  body  rotation.  The  hot  film  was  then  inserted, 
and  the  mean  voltage  vs.  velocity  were  taken,  making  sure  that  the  wake  did 
not  Interfere.  After  three  rotations  the  wake  caused  a downward  drift  In 

voltage.  The  calibration  was  done  for  velocities  from  zero  to  28  cm/sec  and 
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for  salt  water  densities  from  1.000  gm./cc.  to  1.116  which  is  greater  than 
the  range  encountered  in  the  stratified  flow.  The  heat  transfer  was  found 
to  have  no  dependence  on  salinity  as  M.  D'Allonnes,  (1970)  also  found.  It 


was  determined  that  the  heat  transfer  for  a given  film  and  a given  run  was, 
a,  •times  conS+«**'t" 

where  e is  the  voltage  output  from  the  anemometer  bridge  and 
Rw  is  the  resistance  of  the  hot  film. 

If  the  fluid  temperature,  Ta,  and  the  overheat  ratio,  were 
held  constant,  the  heat  transfer  went  as  the  square  root  of  the  velocity. 

The  Rayleigh  No.  of  the  hot  film  was  small  and  the  zero  velocity  transfer  was 
conductive.  These  results  are  shown  in  figure  7.  Immediately  prior  to  each 
stratified  run,  the  hot  film  was  calibrated  in  the  neutral  turbulent  flow  by 
a pitot  tube.  The  ambient  temperature,  Ta,  and  overheat  ratio,  £T,  were 
frequently  checked  and  found  to  be  constant  for  the  results  presented.  Repro- 
ducibility of  velocity  was  generally /n- 0.5% 

e.)  Processing  of  the  Data 

In  processing  the  data  two  criteria  had  to  be  met.  The  electro- 
nics had  to  respond  to  frequencies  from  near  zero  to  5 KHz.  The  spectrum  of 
turbulent  water  occurs  at  lov/er  frequencies  than  that  of  air.  Also  the  system  had  to 
be  on  line,  giving  the  resultant  mean  and  fluctuating  voltages  immediately.  « 
small  drift  in  the  hot  film  voltage  produces  a large  drift  in  the  velocity 
calibration  curve  as: 

e*.  e,*cu.«o)  ^ ^ e*' 

where  is  a constant  times  the  zero  velocity  heat  transfer,  and  B is  the 

slope  of  the  calibration  curve.  This  drift  could  be  picked  up  immediately  by 
observing  the  mean  voltage  output.  Small  air  bubbles  attached  to  the  film 
produce  a large  postive  change  in  he  rms  of  the  fluctuating  voltage.  Because 
of  these  requirements,  a tie  in  with  an  EAI  analog  computer  two  building  away 
was  installed.  A closed  circuit  television  link  was  also  installed  to  observe 
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the  operation  of  the  computer  and,  in  particular,  the  mean  voltage  off  of  a 
four  place  digital  voltmeter  on  the  computer.  An  audio  link  allowed  the 
immediate  observation  of  computer  malfunction  via  the  overload  alarm,  as  well 
as  the  beginning  of  each  integration  period.  A one  minute  integration  time 
was  sufficient  to  include  the  low  frequency  energy  in  the  rms  fluctuation 
output.  As  a result  of  the  visual  and  audio  links,  the  computer  was  made  a 
slave  station  to  the  laboratory. 

In  the  first  eight  months  of  operation  the  Beldon  shielded  cable 
was  adequate  to  avoid  any  noise  pick-up  between  the  laboratory  and  the  computer. 
This  changed  abruptly,  and  the  signal  to  noise  ratio  became  less  than  one. 

After  much  effort  to  locate  the  source,  it  was  decided  to  build  a common 
mode  rejection  device  at  the  computer.  This  eliminated  the  noise  to  less  than 
one  millivolt,  quite  adequate  for  our  purposes.  The  system  v/as  checked  by 
introducing  a known  DC  voltage  and  sine  wave  into  the  transmission  line  at 
the  laboratory  station.  For  frequencies  from  .01  Hz  to  5 KHz  and  amplitudes 
from  .01  volts  to  0.5  volts,  the  computer  output  was  exact  to  within  the 
ability  to  measure  the  input.  This  input  more  than  encompassed  the  range 
found  in  water  born  turbulence.  The  schematic  is  shown  in  figure  2. 

The  mean  velocity  was  determined  by  using  the  calibration  curve. 

That  is; 

a = f 

The  fluctuation  rms  velocites  were  found  by  using  the  slope  of  the  calibra- 
tion curve.  That  is: 

fff  =>  ^ r 

These  profiles  along  with  the  density  profiles  were  then  fed 
into  the  numerical  program  of  the  first  integral  of  the  equations  of  motion, 
which  yield  the  momentum  and  mass  fluxes.  This  is  presented  in  detail  in 
section  four. 


(8) 


3.  The  Neutral  Case 


a.)  Comparison  of  the  Mean  and  Fluctuating  Profiles  with  those  of 
Air  Borne  Turbulence 


In  approaching  the  measurement  of  turbulence  in  a stratified 
flow,  the  investigation  of  water  born  turbulence  in  the  neutral  case  was 
felt  to  be  a necessary  first  step.  This  was  done  to  reproduce  the  well  known 
results  from  air  borne  turbulence,  and  to  verify  the  velocity  instrumentation. 
Ideally,  isotropic,  homogeneous  turbulence,  where  there  is  an  appeal  to 
theory,  would  be  the  first  step.  However,  we  proceeded  directly  to  the  shear 
flow-  the  analog  in  water  of  Laufer's  wind  tunnel  experiment.  This  is  the 
initial  condition  for  the  flow  which  then  becomes  stratified.  The  neutral 
case  has  three  parts:  the  structure  of  a fully  developed,  steady  state,  two 

dimensional  turbulence  shear  flow;  the  requirements  for  two  dimensionality 
and  steady  state;  and  the  nature  of  the  flow  as  the  free  surface  is  approached 
This  region  is  distinct  from  Laufer's  central  region  and  little  is  known  of  it 
The  neutral  case  exhibited  the  expected  "Law  of  the  Wall"  mean 
velocity  distribution.  That  is: 


For  a Re. No.  = 9,400  , figure  S,  k = . 35  -/and  A = 4.0  - .6. 

For  a Re. No.  - 20,000,  figure  9,  k = .41  - .01,  and  A = 5.8  - .8.  Laufer, 

(1951)  for  Re.  No.'^lO^  found  k = .33  and  A=  5.5.  For  large  Re.  No.  10^) 

the  literature  generally  give  k = .41.  There  is  less  agreement  on  the  value 

of  A.  Clauser,  (1956),  found  k = .41  and  A = 4.9  for  Re.  No.  ^ 10^.  All 
of  the  above  results  are  for  hydraulically  smooth  surfaces.  The  logarithmic 
law  of  the  wall,  the  local  equilibrium  region,  extends  over  60%  of  the  depth, 
which  is  approximately  what  Laufer  found,  but  well  beyond  the  theoretical 
prediction  and  what  other  investigators  have  found  in  regards  to  turbulent 
boundary  layers.  The  region  where  the  velocity  gradient  is  linear,  the 
laminar  sublayer,  was  found  to  extend  out  to  where  which  v<as  in 

agreement  with  Laufer.  This  region  is  shown  in  figure  10  and  11.  The  dis- 
tribution of  is  in  agreement  with  Laufer  as  shown  in  figures  12  and 

13.  Initially  a vertical  probe  was  used  causing  a flow  down  the  probe  in 
in  the  region  of  strong  shear  and  causing  the  rms  results  to  be  low  (fig  12); 
the  inclined  probe  shown  in  figure  5 removed  this  problem. 
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b.)  Requirements  for  Two  Dimensionality 


The  mean  velocity  profile  is  quite  insensitive  to  deviations 
from  two  dimensionality  and  steady  state  conditions.  The  most  sensitive 
quantity  to  deviations  from  two  dimensionality  is  the  normalized  rms  value 
of  the  velocity  fluctuations,  For  aspect  ratios  (flow  width  to 

depth)  less  than  4.5,  the  normalized  results  did  not  agree  with  Laufer's. 

For  aspect  ratios  greater  they  did  except  near  the  free  surface.  This  is 
shown  in  figure  12.  The  stream  wise  development  of  the  flow  had  no  effect 
on  the  mean  or  fluctuating  quantities.  For  a ratio  of  inertial  terms  to 
pressure  gradient  of  18%  the  results  were  unchanged. 

c.)  Requirements  for  Steady  Conditions 

The  flow  is  steady  state  when  the  pressure  gradient  balances 
to  Z component  of  the  shear  stress  at  all  points  in  the  flow.  In  a develop- 
ing flow  it  was  found  that  the  effect  of  the  inertial  terms  could  be  accurate- 
ly determined  by  assuming  a logarithmic  velocity  profile  as  shown  in  figure  9 
over  the  entire  cross-section.  The  integrated  streamwise  momentum  equation 
then  becomes: 

where  Fr.  No.  = •fET  , H = depth  of  the  flow.  ^Iso?^  is 

the  slope  of  the  free  surface,  and^j^  is  the  slope  of  the  bottom  of  the 
channel.  TW  is  the  stress  at  the  wall  and  Tf  is  the  stress  at  the 
surface  (about  2%  of  "tW).  The  term  in  bracket  U is  nearly  con- 
stant at  0.78.  Laufer  found  = .0018  for  all  Re.  Nos.  In  this  case 

it  was  .0020.  The  derivation  of  the  above  is  carried  out  in  Appendix  I,  and 
its  validity  is  shown  in  figure  10,  where  the  inertial  terms  and  pressure 
gradient  then  balance  the  stress  at  the  wall.  The  inertial  terms  result  in 
an  18"  imbalance  between  the  pressure  gradient  and  the  shear.  The  above 
derivation  also  indicated  that  if  the  slope  of  the  channel  bottom  is  adjust- 
ed to  equal  the  slope  of  the  surface  the  flow  will  be  steady  state  for  all 
Fr.  Nos.  This  was  the  condition  for  the  results  presented  in  figure  11  where 
the  Fr.  No.  = .96.  The  pressure  gradient  was  of  the  order 

of  so  care  had  to  be  taken  in  its  measurement.  The  pressure  taps 

shown  in  figure  1 were  connected  to  0.5  inch  glass  tubes.  The  difference 
in  height  of  the  water  columns  was  measured  by  a Titian  Tool  Supply  micro 
measuring  scope.  The  columns  weremadeas  large  as  possible  and  still  remain 
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in  the  field  of  the  scope  and  kept  clean.  A wetting  agent  was  used  to  make 
the  meniscus  as  flat  as  possible.  The  scope  read  to  .0001  inch  and  the  zero 
pressure  gradient  reading  was  reproducible  to  - .0004  inch.  The  reproduci- 
bility of  the  pressure  gradient  (the  slope  of  the  free  surface)  was  about  5% 
for  Re.  No.  = 20,000  runs.  The  slope  of  the  bottom  of  the  channel  was  taken 
with  a micro  bubble  gauge.  This  was  constructed  from  two  ten  inch  pieces 
of  machined  aluminum,  hinged  at  one  end  and  with  a micrometer  shaft  at  the 
other  end.  A precision  bubble  gauge  was  placed  in  the  center  of  the  top 
piece  of  aluminum.  In  taking  a reading  of  the  channel  slope  the  bubble  was 
brought  to  the  center  of  the  level  gauge,  and  the  micrometer  reading  was 
compared  to  that  for  zero  slope.  This  technique  was  accurate  to  .0008  inch 
over  a ten  inch  span.  Initially,  a micrometer  point  gauge  was  used  to  back 
up  the  surface  slope  measurement.  All  this  succeeded  in  doing  was  to  measure 
the  evaporation  rate.  An  interesting  side  result  was  that  the  evaporation 
rate  varied  by  a factor  of  20  from  when  the  channel  was  covered  and  the 
water  surface  dirty  to  when  it  was  open  and  the  surface  clean. 

d.)  Mean  and  Turbulent  Properties  Near  a Free  Surface 

i.)  Surface  Stress  and  Tendency  Back  to  Local  Equilibrium 

The  results  near  the  free  surface  differed  from  those  of  Laufer's 
near  the  center  of  the  wind  tunnel.  Both  the  mean  and  fluctuating  velocities 
increased  beyond  the  Laufer  results  as  the  free  surface  was  approached,  as 
shown  in  figures  14  and  15.  The  significant  difference  between  a free  sur- 
face and  the  center  of  a wind  tunnel  is  the  appearance  of  a molecular  film 
and  a resultant  stress  at  the  free  surface.  This  non  zero  boundary  condition 
is  usually  bothersome  to  the  laboratory  experimentalist,  but  is  of  great 
interest  to  the  oceanographer.  The  wind  exerts  a surface  stress  on  the  ocean 
currents.  For  a zero  surface  stress  flow  the  region  near  the  surface 
(0-t  f ^ ilo  ) is  difficult  to  characterize  because  it  is 
not  in  local  equilibrium.  Diffusion  is  important  and  the  properties  of  this 
region  depend  upon  what  is  happening  elsewhere  in  the  flow.  As  a surface  stress 
is  applied,  the  flow  tends  towards  a constant  stress  Couette  flow.  The  boundary 
condition  downstream  at  the  dam  is  that  the  Froude  Number  be  one.  This  implies 
that  the  flow  remains  at  a nearly  constant  depth.  The  velocity  at  the  surface 
increases  over  the  zero  stress  condition,  and  for  a fixed  flowrate,  conservation 
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of  mass  requires  that  the  velocity  and  shear  near  the  wall  decrease.  The  wall  stress 
decreases  as  the  surface  stress  increases  and  the  flow  tends  towards  local  equili- 
brium where  the  diffusion  of  energy  is  unimportant  and  a balance  between  local 
energy  production  and  dissipation  occurs. 


7) 


ii)  Discussion  of  a Laminar  Sublayer 

The  stress  at  the  free  surface  was  determined  thusly: 

Following  the  same  reasoning  used  at  the  bottom  of  the  channel, 
the  laminar  sublayer  is: 


giving  ^ 

The  equations  of  motion  in  this  region  require  that: 

7 = “ii  - f 


8) 

9) 

10) 
11) 


and 


thus 


and 


r fftls-Cti-ti)* h)* 

Us  (t^*rs)t^[K,M  - (rw-Ts)»ygjk 

These  results  are  plotted  in  figure  14  and  clearly  do  not  conform 
to  the  measured  profile.  This  region  is  turbulent.  The  rub  is  that  tUs  • 
does  not  characterize  the  velocity  here.  If  the  surface  velocity  is  used, 
then  for  iu»/v  , f 

the  laminar  sublayer  occurs  at  ^ ^ 

This  is  considerably  closer  to  the  surface  that  it  was  able  to  measure  and  is 
smaller  than  surface  disturbances.  The  free  surface  is  hydraulically  rough. 
The  surface  stress  is  larger  than  equation 7.)  but  has  an  upper  limit  from 
the  imbalance  of  the  correctei^^^H  and  Tw  . which  is  about  5%  of  1T«#» 
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lii)  Theoretical  Attempts  to  Describe  the  Velocity  Structure 

If  we  assume  that  the  region  near  the  surface  is  turbulent 
and  in  local  equilibrium,  similarity  arguments  can  be  used  to  predict  the 
mean  velocity  structure.  Prandtl,  (1949),  and  Townsend,  (1961)  describe 
an  equilibrium  layer  as  a region  where  there  is  a balance  between  local 
turbulent  energy  production  and  dissipation,  and  where  the  macroscale  of  the 
turbulence  is  small  compared  with  the  scale  of  the  change  in  stress.  Under 
these  conditions  the  formulation  for  equilibrium  layers  is: 

where  Von  Karman's  constant 

.HI 

For  our  co-ordinate  system  and  for  the  free  surface  region  this  becomes: 


where  H = depth  of  the  flow 
Integrating  from  ik  to  h-€we  find: 

'2)  W.{ljr  UCn-S)  - -Xl  ^i-Vh  - ■&'] 


■La  I ■ UI-VI. . -1^ 

In  this  equation  S is  the  point  of  velocity  measurement  closest  to  the  surface. 
The  results  from  this  equation  are  plotted  in  figures  14  and  15  for  qpi6 
and = Figure  14  shows  that  the  predicted  velocity  gradient  is  steep- 


er than  the  measured.  Because  the  diffusion  of  turbulent  energy  is  probably 
still  important  for  these  low  values  of  ^/Tw  y Kwt  is  underestimated. 


Figure  15  indicates  a much  better  agreement  between  theory  and  experiment  for 
this  higher  ration  of  Vr^  The  agreement  is  good  out  tc||I*i^H  = .08 
where  the  theoretical  gradient  again  steepens  beyond  the  measured  value. 


4.  Investigation  of  the  Stratified  Case 

a.)  Equations  of  Motion  and  Mass  Transfer 

The  aim  of  this  portion  of  the  experiment  is  to  determine  the  mass 
and  momentum  fluxes  across  the  channel  for  various  bulk  Ri.  Nos.  and  at  an 
approximately  constant  Re.  No.  Profiles  of  mean  density,  mean  velocity,  and 
the  time  average  of  the  fluctuating  velocity  were  taken  at  cross  sections  one 
foot  apart  in  the  center  of  the  channel  where  y*0  and  ^ =0  from  symnetry. 
This  information  along  with  the  boundary  condition,  the  wall  shear  and  the 
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pressure  gradient,  was  used  in  conjunction  with  the  first  integral  of  the 
equation  of  motion  to  determine  the  fluxes. 

The  equations  of  motion  for  an  incompressible  turbulent  flow  are; 

The  boussinesq  approximation  has  been  used  as 

O OH  for  all  runs, 

where  is  the  difference  in  mean  density  across  the  flow. 

Using  the  normal  procedure,  for  a time  independent  mean  flow,  of  separat- 
ing the  dependentvariables  into  a mean  and  fluctuating  component,  let; 

UL^r  U.4>U*  , r f4>  f' 

Vs/^VW+W'  , P*rP«*P» 

By  definition  the  integral  of  the  fluctuating  quantities  over  a period  of  time 
is  zero.  This  time  period,  t,  is  found  by  experiment.  Thus,  after  taking  the 
integral  over  t,  the  equations  of  motion  become: 


17.) 


s O 


18.) 

19.)  ^ 


<^^^^is  neglected,  as  thedata  from  the  velocity  profiles  of  figure  16  reveals: 


20.) 


Because  the  stream  wise  scale  is  much  greater  than  the  vertical  scale,  that  is 
from  figure  16: 

the  hydrostatic  form  of  equation  (20)  can  be  used. 

Me  then  have: 

The  conservaton  of  salinity  is: 

'X&K*  1s  neglected  as  the  mean  density  profiles  from  figure  17  reveal 
Using  equation  (18)  and  letting  . ,t  ^ . s- 
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23. ) 

24. ) 


The  last  term  in  equation  (24)  cannot  be  measured  with  existing 
toclinii]uos.  However,  if  we  assume  that  the  turbulent  scales,  the  eddies, 
effecting  the  momentum  transfer,  are  representative  of  thoseeffecting  the 
mass  transfer  and  these  scales  go  as  in  the  streamwise  direction  and  ia'W ’ 

in  the  vertical , we  have  from  figures  18  and  19  for  =0.8 

and  ^ “7557?^  ' ' * 

VI  T’ICl 

rFSFvIm  T" 

Thus,  equation  (24)  becomes: 


At.  ^ . •»>V  JFWT 

25)  -sr  ♦'^5$.+ '"■55*  ^Ti«'  -f* 


Next,  we  integrate  equation  (21)  from  zero  to  some  value  of  Z and  take  the 
stream  wise  (x)  gradient.  This  yields  for  equation  (23) 

Where  Pw  is  the  pressure  at  the  wall.  Integrating  equation  (26)  from  zero  to 
Z and  using  equation  (17)  we  have 
27) 


lu  \ we  nave; 

Where  is  the  stress.at  the  wall.  ^ 

Now  let 

’ s Z - Jl'  1 ^ 

* -f  fi-  Ctt^w^54i'  - *' 


Equation  (27)  becomes: 

28)  ra^ 


Equation  (28)  gives  the  distribution  of  the  momentum  flux,  , as  a 

function  of  the  boundary  conditions  and  the  measured  stream-wise  gradients. 
Term  1 on  the  right  hand  side  is  the  stress  at  the  wall;  Term  2 is  the  stream 
wise  pressure  gradient  at  the  wall;  term  3 and  4 are  the  contribution  to  the 
pressure  gradient  by  the  internal  density  gradient  which  will  be  referred  to 
as  the  baroclinic  pressure  gradient.  Term  5 is  the  integrated  divergence  of 
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the  horizontal  momentum  flux  and  term  6 is  the  mean  vertical  momentum  flux. 
The  contribution  of  these  terms  to  for  a Re.  No.  of  10,000  and  a Ri.No. 
of  0.33  is  presented  in  Table  I. 

Now  integrating  equation  (25)  and  using  equation 

(17): 


29. 


• • • 


As  the  flow  is  recirculating,  there  is  a slow  increase  in  the  ambient  density 
in  time.  However,  in  the  region  where  the  salt  flux  is  important,  that  is 
for  */K  , figures  20  and  21  reveal: 


Also,  as  is  shown  in  Table  II; 


Thus,  equation  (29)  becomes: 


w/ 
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Equation  (30)  gives  the  distribution  of  the  turbulent  salt  flux, 
across  the  flow  as  a function  of  the  measured  stream  wise  gradients.  Term  1 on 
the  right  hand  side  is  the  integrated  mean  h.rizontal  flux  divergence  and  2 is 
the  mean  vertical  flux.  The  contribution  of  these  terms  to  for  a Re. 

No.  of  10,000  and  a Ri.  No.  of  0.33  is  presented  in  Table  I. 

If  we  let: 


and  transform  equation  (28)  and  (30)  to  numerical  form,  we  then  have: 


31)  T(«)s  CVtK.a  ; Tw«  V 

.9!: 

AH-  ’ «,=l  A1«.  • 

The  fluxes  were  then  determined  by  using  the  raw  data  from  the  measured  profiles 
and  programing  equations  (31)  and  (32)  for  a digital  computer.  The  Richardson 
Number:  _ 

was  determined  in  the  same  manner.  The  program  also  yielded  Ks  and  Km  by 
dividing  the  fluxes  by  the  measured  vertical  gradients. 

Because  of  the  scatter  in  as  shown  in  figure  18,  its  stream 
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wise  gradient  was  neglected.  However,  as  figures  16  and  18  indicate  for;£  /H=0.17 

It  was  generally  much  smaller  than  this. 

Also,  the  scatter  in  measuring  the  pressure  gradient  was  ^ 50%  (about 
five  times  poorer  than  in  the  neutral  case).  Consequently,  the  pressure  gradient 
was  adjusted  to  provide  zero  stress  at  the  free  surface. 

The  pressure  gradient  was  a significant  term  in  the  momentum  flux  equation  as 
indicated  in  Table  I.  As  a result,  large  scatter  in  the  pressure  gradient 
created  large  scatter  in  the  resultant  momentum  flux. 

b.)  The  Measured  Profiles 

i.)  Gross  Features  of  the  Developing  Flow 

Runs  were  made  at  a Re. No.  of  ^ 10,000  and  Ri.Nos.  of  0.22,  0.30, 

0.33,  0.36,  and  0.38.  The  Re.  No.  was  defined  by  using  the  maximum  velocity 
and  the  depth  of  the  flow.  The  bulk  Ri.  No.  was  defined  by  using  the  maximum 
velocity,  the  depth  of  the  flow,  and  the  density  difference  between  the  top 
and  bottom  of  the  flow. 

The  earlier  runs  made  at  Ri.  Nos.  of  0.22,  and  0.33  had  considerable 
scatter  in  regards  to  the  density  profiles.  Consequently,  the  data  was  not 
reduced  and  the  bulk  Ri.  Nos.  must  be  considered  approximate.  The  mean  and 
fluctuating  velocity  profiles  are  useful  and  are  presented  in  Appendix  II. 

The  measured  pre'iles  presented  in  the  graphs  at  the  end  of  this  sec- 
tion and  in  Appendix  II  indicated  that  for  a given  run  over  most  of  the  cross 
section  the  mean  flow  is  accelerating  and  U/*  , the  wall  stress  and  the 
bulk  Ri.  No.  are  decreasing  in  IP  . If  the  channel  were  long  enough,  the  flow 
would  return  to  the  neutral  case  and  the  classical  turbulent  flow  presented 
in  figures  22  to  23.  At  some  point,  the  possible  critical  Ri.  No.,  the  stream- 
wise  gradients  of  uJ*  and  must  become  positive  to  return  to  the  neutral 

conditions. 

In  going  from  one  run  to  another  of  lower  bulk  Ri.  No.  the  opposite 
occurs.  The  wall  stress  increases, the  maximum  velocity  decreases.  The  flow 
is  tending  baci  towards  neutral  conditions. 

ii)  The  mean  Velocity  Results 

The  mean  velocity  results  were  of  high  quality,  reproducibility 
in  leneral  being  better  than  one  percent  as  indicated  in  figures  16,  II-l  and 
II  - 3.  More  than  twice  as  many  runs  were  made  than  are  presented  here.  Ini- 
tially a pitot  tube  was  used  to  make  a survey  over  a wide  range  of  Ri.  Nos.  to 
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determine  what  ranqe  *^se  the  more  detailed  hot  film  instrumentation.  The 
pitot  tube  proved  to  be  useless  in  the  stratified  case  as  the  scatter  shown  in 
figure  6 indicates.  The  hot  films  were  not  of  uniform  quality  and  had  a shelf 
life,  particularly  after  having  been  used  once  in  saline  water.  Sometimes  they 
would  deteriorate  abruptly  as  was  the  case  in  the  third  cross  section  of 
figure  II-l.  Other  times  they  would  exhibit  a slow  drift.  As  the  mean  velo- 
city goes  as  the  fourth  power  of  the  measured  voltage,  an  on-line  system  with 
a very  accurate  voltmeter  is  necessary  to  continuously  monitor  the  equipment 
performance.  A closed  circuit  television  system  read  the  voltage  output 
directly  off  of  the  four  place  digital  voltmeter  on  the  computer.  When  the 
hot  films  worked,  they  worked  quite  well. 

iii)  The  fluctuating  Velocity  Results 

The  rms  of  the  voltage  fluctuations  was  also  continuously  moni- 
tored. The  results  of  the  run  shown  in  figure  1 1-2  are  presented  even  though 
a significant  change  in  the  calibration  of  the  hot  film  occurred.  The  rms 
fluctuations  depend  upon  the  slope  of  the  calibration  curve  and  are  less  affect- 
ed by  a displacement  of  it. 

Another  reason  for  the  on-line  system  was  to  immediately  detect 
and  reject  the  rms  results  that  were  affected  by  air  bubbles  forming  on  the 
hot  film.  The  rms  of  the  fluctuating  voltage  would  in  general  show  an  immedi- 
ate change  when  this  occurred.  Air  bubbles  coming  out  of  solution  and  adhering 
to  the  back  of  the  hot  film  caused  virtually  no  problem  in  the  neutral  fresh 
water  runs  when  the  water  was  sufficiently  aged.  However,  in  the  stratified 
case,  particularly  in  the  interior  regions  where  strong  mixing  occurred,  this 
was  a problem.  The  dissolving  of  salt  into  fresh  water  lowers  its  saturation 
level  and  makes  air  bubbles  more  prone  to  come  out  of  solution.  Small  bubbles 
would  reduce  the  mean  heat  transfer,  a large  bubble  would  increase  it,  and 
sometimes  two  large  bubbles  close  to  each  other  would  reduce  it.  All  bubbles 
seemed  to  increase  the  fluctuating  component.  The  scatter  in  the  rms  of  the 
velocity  fluctuations  seened  to  come  from  bursts  of  energy.  Intermitancy  is  a 
problem  found  by  many  other  investigators.  There  is  also  the  possibility  that 
small  bubbles  , formed  towards  the  end  of  the  integration  period,  went  undetected 
and  had  some  effect.  The  mean  conditions  remained  constant  as  monitored  by  the 
flow  meter,  the  depth  and  temperature  of  the  flow,  and  the  reproducibility  of 
the  mean  velocity.  Also  the  salt  water  injection  rate  proved  to  be  quite  constant. 


An  interesting  qualitative  observation  in  regards  to  the 
fluctuations,  is  that  in  the  region  where  they  are  largest,  internal  waves 
were  visually  observed. 

iv)  The  Density  Profiles 

Considerably  more  difficulty  developed  in  measuring  the  mean 
density  profiles.  The  problems  were  not  wholly  resolved  until  the  last  run, 
represented  in  figures  16  to  33;  consequently  the  previous  runs  should  be 
viewed  qual itatively.  Because  the  flow  is  recirculating,  there  is  an  increase 
in  the  ambient  density  of  the  flow  upstream  of  the  salt  injec- 

tors. This  rate  of  increase  is  negligible  in  the  equations  of  motion  as  was 
previously  shown,  but  must  be  taken  into  account  in  correcting  the  measured 
density  profiles.  The  mean  conditions  of  the  flow  remained  constant,  as  did 
the  injection  flow  rate  of  salt  water  shown  in  figure  3;  however,  the  rate  of 
increase  of  the  ambient  density  exhibited  considerable  scatter  until  the  last 
run.  Initially,  the  ambient  density  was  monitored  by  placing  the  density  probe 
near  the  surface  and  downstream  of  the  test  section  during  the  three  hours  the 
velocity  profiles  were  taken.  This  was  found  to  indicate  a much  higher  rate 
of  ambient  density  increase  than  was  true.  The  wake  of  the  velocity  probe  caused 
this  anomaly.  In  the  last  run  the  density  near  the  surface  was  monitored  in  the 
test  section  when  the  density  profiles  were  taken.  The  results  are  shown  clearly 
in  figure  20. 

Another  possible  source  of  error  that  was  corrected  in  the  last 
run,  figures  16  to  33,  was  the  possible  mixing  of  fluid  from  one  data  point  when 
the  next  data  point  was  taken.  Rather  than  use  the  hydrometer  as  a continuously 
reading  instrument  as  i t was  designed  for,  it  was  emptied  after  each  data  point 
and  slowly  filled  at  the  next  one. 

In  using  the  density  probe  ideally  the  velocity  of  fluid  entering 
the  probe  should  match  or  be  less  than  that  of  the  flow  at  that  point.  In  the 
last  run  the  maximum  possible  filling  time  of  the  hydrometer  of  four  minutes 
was  used.  This  limit  is  imposed  because  of  the  storage  capacity  of  saturated 
salt  water  and  the  time  required  to  take  all  of  the  data  which  was  about  six 
hours.  The  probe  entry  velocity  for  a four  minute  filling  time  was  42  cm/sec. 

In  an  unstratified  flow  the  probe's  resolution  at  an  approach  velocity  of  18  cm/sec 
was  observed  by  the  use  of  dye  streaks.  The  resolution  was  essentially  that  of 
the  probe's  dimension  of  .089  cm.  This  technique  was  not  able  to  be  used  in  the 
low  velocity,  strong  gradient  region  near  the  wall  in  the  stratified  case, 
because  of  the  fluctuating  density  and  index  of  refraction. 
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This  is  a source  of  possible  systematic  error  in  the  density  profiles.  It 
would  affect  the  vertical  gradients  near  the  wall,  but  should  have  little 
effect  on  the  streamwise  gradients  as  the  velocity,  the  density  and  their  ver- 
tical gradients  are  similar  at  each  cross  section. 

c.)  The  Results  of  the  First  Integral  of  the  Equations  of  Motion 

i.)  The  Momentum  and  Mass  Fluxes 

Initially,  to  determine  the  momentum  and  mass  fluxes,  three  cross 
sections  were  taken.  These  mean  and  fluctuating  profiles  are  shown  in  Appendix 
!’  . This  was  reduced  to  a two  cross  section  run  so  that  more  time  would  be 
allowed  for  individual  measurements,  particularly  in  regard  to  the  density  pro- 
files. The  numerical  program  for  equations  (31)^  (32)  was  worked  out  manually 
in  its  entirety  for  six  of  the  16  values  of  n in  the  case  of  the  mass  transfer, 
and  for  two  values  of  n in  the  momentum  case.  Also  several  of  the  individual 
terms  in  the  momentum  equation  for  larger  values  of  n were  checked,  and  the  first 
term  in  the  mass  transfer  equation  was  worked  out  for  all  values  of  n. 

The  error  bars  were  determined  by  introducting  in  the  computer 
program  an  error  of  2 1/2"  in  the  velocity  (an  error  of  12%  in  the  streamwise 

gradient)  at  the  lowest  value  of  n where  large  changes  occur;  that  is  for  n=4 

and  Z/h=.062.  This  error  is  considerably  larger  than  the  measurement  errors 
at  the  two  adjacent  points  that  determine  the  vertical  gradient.  However,  be- 
cause the  velocity  used  to  determine  the  streamwise  gradient  had  to  be  read 
from  the  curve  drawn  connecting  the  two  adjacent  points,  the  uncertainty  in 
curvature  at  this  value  of  n might  produce  this  magnitude  of  error.  Similar 
errors  in  velocity  were  introduced  at  larger  values  of  n.  They  produced  much 
smaller  error  bars.  E>"rors  in  density  are  not  included  as  they  also  have  a 
much  smaller  effect.  The  inertial  terms  dominate  as  Table  I indicates  and  the 
percentage  error  in  the  streamwise  density  gradient  is  smaller.  The  error  bars 
in  the  mass  transport,^*  , are  larger  as  two  relatively  large  terns  that  are 

nearly  equal  are  subtracted  from  one  another.  This  is  pointed  out  in  Table  1 

and  r i qure  24 . 

A strong  indicator  of  the  reliability  of  the  salinity  flux  results 
is  the  conservation  of  mass  at  the  two  cross  sections,  that  is: 
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The  difference  between  these  is  the  first  term  in  the  mass  flux  equation  and 
should  be  nearly  zero  at  the  surface.  In  the  last  run,  figures  16  to  33,  these 
terms  were  within  0.08^,  of  each  other.  In  the  previous  runs  they  were  from  1% 
to  4T  of  each  other.  The  preceeding  error  analysis  is  not  intended  to  give 
the  absolute  error  bars,  but  rather  to  give  an  indication  of  what  they  are. 


ii)  Energy  Budget  for  a Developing  Flow 

Considering  the  momentum  flux,  and  particularly  that  component 
produced  by  turbulence,  Wi'W*  , of  figures  19  and  26,  the  most  striking  fea- 
tures is  that  it  is  negative  over  the  lower  portion  of  the  flow.  At  this 
point,  it  is  useful  to  remember  that  in  a turbulent  system  we  have  a balance 
between  the  mean  advection  and  pressure  gradient  and  the  fluctuating  iner- 
tial tenns  which  sometimes  act  as  dissipative  terms.  It  is  clear  from  the 
results  of  ul*  and  that  if  we  define 

^ --ST. , 

then,  for  the  upper  portion  of  the  flow  where 
then, 

O 

In  regard  to  u*w*  the  turbulent  energy  budget  equation  should  be  reconsidered 
in  light  of  the  developing  aspect  of  the  flow.  In  general  from  Phillips,  (1966), 
this  is: 


(33)  1 0^  ^ ^ 
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Where  4 is  viscous  dissipation  and  is  positive.  For  a two  dimensional  devel- 
oping flow  this  becomes: 

(34)  , ^ y w* - 4 > 
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if  we  let: 


we  then  have: 
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In  the  steady  state  case  the  right  hand  side  of  equation  (35)  is  zero,  limit- 
ing the  Rf.  No.  to  less  than  one.  In  our  case  the  first  two  terms  on  the 
right  represent  the  advection  of  turbulent  energy.  The  third  and  fourth 
terms  represent  the  vertical  and  streamwise  diffusion  of  turbulent  energy. 

The  pnysi.dl  meaning  of  the  remaining  terms  is  not  clear.  Thus,  in  a devel- 
oping flow  the  Rf.No.  can  be  greater  than  one  and  turbulence  still  persist 
if  the  vertical  and  streamwise  gradients  of  the  turbulent  fluctutations  are 
negative.  To  say  this  another  way  the  Rf.  No.,  as  defined,  does  not  properly 
assess  the  kinetic  energy  available  for  mixing.  Rearranging  equation  (35) 
we  have; 
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In  equation  (36)  the  Reynolds  stress,  the  term  on  the  left,  is  balanced  by 
the  bouyancy  and  viscous  forces,  the  first  two  terms  on  the  right,  as  well 
as  the  intertial  forces  of  turbulent  advection  and  diffusion,  the  third, 
fourth,  fifth  and  sixth  tenns  on  the  right.  The  physical  meaning  of  the  last 
term  is  not  clear.  Thus,  and  V(w%^  (as  is  always 

positive)  can  be  negative  if  the  vertical  and  streamwise  gradients  of  the  tur 
bulent  fluctuations  are  strongly  negative. 

The  implications  of  equation  (36)  can  be  seen  in  the  results  pre 
sented  in  figures  19  and  27.  The  Reynolds  stress,  - . and  the  Aus- 

tausch  coefficient,  , are  negative  in  the  region  near  the  wall.  One  of 
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the  terms  on  the  right  hand  side  of  equation  (22)  can  be  approximately  deter- 
mined from  the  measured  profiles.  For  */h  =.20  of  figures  18  and  19  we  have: 

lA,' (-2  0)  s.  —.066  tt-'/jcc*- 

_.o> 

The  downstream  station  four  results  should  be  considered  more  reliable  as  the 
streamwise  gradients  used  to  determine  them  more  properly  apply  here.  The  impli- 
cations of  equation  (35)  can  be  seen  in  the  upper  region  of  the  flow.  The  flow 

is  highly  turbulent  (Km^  is  ^ 15  times  that  of  the  viscosity  as  shown  in  fig- 
ure 27)  and  the  Ri.  No.  is  of  order  three.  The  Rf.  No.  involves  the  ratio  of 
ks/kmz  and  is  possibly  smaller  than  the  Ri.  No.  Unfortunately,  Ks  is  not  pre- 
sented for  this  region  as  both  the  flux  and  the  gradients  approach  zero,  making 

Ks  indeterminate.  However,  because  this  region  is  quite  turbulent,  Ks  would  be 
similar  to  Kmz  and  the  Rf.  No.  would  be  close  to  the  value  of  the  Ri.  No. 

iii. )  Austausch  Coefficients  and  the  Local  Richardson  Number 

Figure  28  shows  the  momentum  and  salinity  coefficients  as  a 
function  of  the  local  Ri.  No.  for  the  upstream  and  downstream  cross  sections 
of  the  last  run.  The  momentum  coefficient  shows  consideraMn  scatter  and  does 
not  show  the  expected  reciprocal  relationship  to  the  local  Ri.  No.  as  it  is 
defined.  It  is  not  a single  valued  function  of  the  local  Ri.  No.  This  is 
thought  to  be  the  result  of  the  importance  of  advection  in  the  momentum  bal- 
ance in  the  different  regions  of  the  flow.  This  will  be  discussed  in  more 
detail  in  the  next  section.  The  salinity  coefficient  behaves  in  the  expected 
manner  with  the  coefficient  decreasing  as  the  local  Ri.  No.  increases  with 
the  exception  of  the  region  near  the  wall.  In  going  from  near  zero  to  a 
Ri.  No.  of  3.25  the  coefficient  decreases  by  a factor  of  five,  although  it  is 
still  some  200  times  its  molecular  value.  For  a Ri.  No.  of  3.25  the  slope  of 
the  curve  is  still  significant.  This  indicates  that  a much  larger  local  Ri.  No. 
is  necessary  to  determine  a possible  asymototic  value. 

iv. )  Regimes  of  the  Flow 

By  considering  the  various  terms  that  contribute  to  the  momentum 
flux  as  a function  of  depth  three  regimes  of  the  flow  can  be  established.  In 
figure  29  the  mean  advective  fluxes  and  their  sum  are  plotted.  Taken  individu- 
ally they  are  mathematical  book  keeping.  It  is  their  sum  which  has  physical 
meaning  and  is  plotted  In  figure  30  to  determine  the  regimes  of  the  flow.  For 
z/h<.35  the  Reynolds  stress  is  negative  and  is  balanced  by  the  sum  of  the  mean 
advective  fluxes.  The  momentum  Austausch  coefficient  is  negative  and  of  the 
order  of  the  molecular  viscosity  as  shown  in  figure  31.  Also  from  figure  16 

(23) 


the  streaniwise  mean  velocity  yradient  is  negative.  Figure  18  shows  that  the 
streamwise  fluctuating  velocity,  Vfc'*  , has  a greater  negative  gradient  here 
than  in  tlie  rest  of  the  How.  This,  considering  the  energy  budget  equations 
is  consistent  with  being  negative.  At  z/H  = .35,  the  curves  of  W *■ 

at  station  two  and  station  four  intersect  and  the  gradient  is  zero.  Another 
feature  of  tin’s  region  is  that  the  vertical  mean  velocity  gradient  with  the 
exception  of  very  close  to  the  wall  is  increasing  in  a stream  wise  direction. 
This  is  consistent  with  u'  and  w'  being  greater  than  zero  or  both  being  less 
than  zero.  Figure  16  shows  that  the  mean  velocity  profiles  intersect  at  z/H 
^ .30  and  above  this  point  the  flow  is  accelerating.  This  along  with  the  fact 
that  internal  gravity  waves  are  observed  in  the  lower  regime  of  the  flow  sug- 
gests that  the  critical  layer  mechanism  is  transferring  momentum  from  the 
fluctuations  to  the  mean  flow.  Bretherton  (1967,  1969j[  points  out  that,  if 
the  local  Ri.  No.  is  greater  than  0.25,  internal  gravity  waves  are  attenuated 
by  a factor  exp  Klo - and  their  momentum  is  transferred  to 

the  mean  flow.  Figures  32  and  33  indicate  that  the  local  Ri.  No.  is  approxi- 
mately 0.25  at  the  point  where  the  flow  begins  accelerating.  For  .35<z/H<.85 
the  Reynold  stress  is  large  and  positive  and  is  balanced  by  the  baroclinic 
pressure  gradient.  That  is; 

t §* 

-■TV. 

The  momentum  Austausch  coefficient  is  large  and  positive  as  is  the  local  Ri.No. 
For  the  surface  region,  Z/h  > .85,  the  Reynolds  stress,  the  Austausch  coeffic- 
ient, and  the  local  Ri.  No.  all  tend  to  zero. 

The  salinity  flux  across  the  flow  exhibits  three  regimes  also, 
although  they  are  not  coincidental  with  the  momentum  case.  Figure  24  shows 
the  mean  advective  fluxes  that  contribute  to  the  turbulent  salinity  flux.  They 
are  both  large  compared  to  the  turbulent  component  and  are  of  opposite  sign. 

For  the  wall  region  Z/h <.06,  Figure  28  shows  that  the  salinity  Austausch 
coefficient  is  small  and  increasing,  although  two  orders  of  magnitude  greater 
than  the  molecular  value.  The  Ri . No.  is  also  small  and  increasing.  For  the 
central  Region,  .06  < Z;':<0.6,  the  Austausch  coefficient  is  large  and  exhibits 
the  expected  reciprocal  relationship  to  the  local  Ri.  No.  Also  shown  in  Figure 
28  is  the  empirical  relationship  where  is  the  value 

for  near  neutral  conditions.  This  relationship  has  been  used  with  r=8  by 
several  investigators  to  model  the  thermocline  and  shows  a rather  poor  fit  to 
the  data.  A curve  wither  adjusted  to  six  to  fit  the  data  at  Ri.  No.  = 3.0 
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does  not  fare  much  better.  The  slope  of  the  data  at  Ri.  = 3.25  is  still 
significant  and  is  larger  than  the  two  above  mentioned  empirical  recipes. 

For  this  reason  is  did  not  seem  worthwhile  to  develop  an  empirical  formula 
for  kj  at  large  Ri.  Nos.  And  finally,  for  Z/h  > .6,  the  flux  and  the 
density  gradient  go  to  zero,  making  the  Austausch  coefficient  indeterminate. 
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FIG.  i3  NORMALIZED  RMS  FLUCTUATING  VELOCITY  vs.  NORMALIZED  DEPTH 
NEAR  THE  WALL  FOR  DIFFERENT  ASPECT  RATIOS  AND  Re.  Nos. 
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FIG.  21  MOMENTUM  AND  SALINITY  FLUXES  vs.  NORMALIZED  DEPTH 
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FIG.  26  REYNOLDS  STRESS  vs.  DEPTH  UPSTREAM  STATION  2.  Re.  No.  = tO.OCC  ; Ri.  No.  *.38 
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FIG.  2 7 MOMENTUM  AND  SALINITY  AUSTAUSCH  COEFFICIENTS  vs.  DEPTH 


MOMENTUM  AUSTAUSCH  COEFFICIENT  vs.  Ri.  No.  SALINITY  AUSTAUSCH  COEFFICIENT  vs. 


-Ol 


04 


FIG.  32  LOCAL  Ri  No.  vs.  NORMALIZED  DEPTH  DOWNSTREAM  STATION  4 
Re.  No.  = 10,000  ; Ri.  No.  ^ .33 
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Contributions  to  the  Momentum  and  Salinity  Fluxes  from  the 
Various  Terms  in  the  Eouatlons  of  Motion 
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APPENDIX  I 


Requirements  for  a Steady  State  Flow  In  the  Neutral 
Case  Using  a Logarithmic  Velocity  Profile 


The  geometry  of  the  flow  Isi 


From  continuity  we  havei 
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where  Z Is  the  point  of  the  mean  velocity  between  H, 
and  at  x,  + 
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combining  (2)  and  (3)  we  have. 
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To  first  order  we  then  have. 


The  momentum  equations  are j 
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Integrating  equation  (7), 
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equations  (6)  and  (9)  yield 
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Thus  equation  (10)  uecomes 

00  U ^ ^ ^ w . <1  c csC  ( iU  -t.,,  G)  r Vi  Tf>,  Z) 


.teJ 


Using  the  continuity  equation i 
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and  Integrating  from  0 to  H,  equation  (12)  becomes 
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Now  as  U)  = ^ <<  1 


We  have 
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Letting, 
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we  have  I 
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Using  equation  (5)  -^U»  /I*)  d 2: 


we  have 
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Where  the  value  of  the  brackets 
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Is  nearly  constant 


at  0,78.  Laufer,  ref  15,  found  the  ratio  of  ~ .0018 

U» 

for  all  Re. Nos.  The  value  here  Is  .0020 


and 


Thus 


The 


Thus 


Thus 


Using  the  data  from  figure  10 

p p 
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corrected  pressure  gradient  balances 
Using  the  data  from  figure  H > 
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, by  tilting  the  channel  the  flow  is  steady  state. 


APPENDIX  II 


Preliminary  Mean  and  Fluctuating 
Velocity  Profiles  at  Various 
Richardson  Numbers 
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FIG.  II-l 


1 1- 3 MEAN  VELOCITY  vs.  DEPTH 


FIG.  1 1 -4  RMS  FLUCTUATING  VELOCITY  vt.  DEPTH 


